Interaction effect on topological classification of superconductors in two dimensions 
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We introduce a new class of superconductors (SCs) in two spatial dimensions with time reversal 
symmetry and reflection (i.e., mirror) symmetry. In the absence of interactions, topological classes 
of these SCs are distinguished by an integer-valued (Z) topological invariant. When interactions are 
included, we show that the topological classification is modified to Zg. This clearly demonstrates 
that interactions can have qualitative effect on topological classifications of gapped states of matter 
in more than one dimension. 



Introduction: Topological phases are fully quantum 
mechanical states of matter which are not characterized 
by classical symmetry breaking. While gapped in the 
bulk, quite often, they are accompanied by gapless ex- 
citations at their boundary, signaling highly entangled 
nature of their ground state. Since the discovery of the 
integer quantum Hall effect (IQHE) [Q, the list of topo- 
logical phases in nature has been expanded, in particular, 
by the recent discovery of topological insulators in two 
and three dimensions (2D and 3D) in systems with strong 
spin-orbit coupling P-[lH, and the identification of '^He 
B as a topological SC (superfluid) ||l2|. Unlike the IQHE, 
the topological character of these topological insulators 
and SCs (i.e., the stable gapless edge or surface modes) 
is protected by time-reversal symmetry (TRS). The pres- 
ence or absence of a topological distinction among gapped 
phases for a given set of symmetries and for given spatial 
dimensions can be studied systematically, and is summa- 
rized in the "periodic table" of topological insulators and 
SCs for non-interacting fermions [p"4|-p^. 

Since interactions are ubiquitous in real materials, 
a natural question is whether and how interactions 
could modify topological classifications obtained for non- 
interacting systems. In other words, are there cases 
where interactions qualitatively affect topological clas- 
sifications? Even though this question has not been fully 
or thoroughly answered, an exciting progress has been 
made in Ref. in which Fidkowski and Kitaev ex- 

plicitly showed that the putative Z classification of one- 
dimensional (ID) "non-interacting" SCs with unusual 
time reversal symmetry ("BDI class") is modified to Zg 
when interactions are included. This Zg classification of 
the ID BDI SCs with interactions is further illustrated 
from the study of entanglement spectrum ||l^, |l^. In 
more than one dimension, stability of topological classifi- 
cations in a number of classes of non-interacting topolog- 
ical insulators and SCs is argued from the perspective of 
low-energy topological response theory jl^, . For 

bosonic (or spin) systems, a systematic approach of con- 
structing symmetry protected topological phases in gen- 
eral spatial dimensions was recently proposed [|2j . How- 
ever, whether and how interactions in fermionic systems 



in 2D and 3D can dramatically modify their topologi- 
cal classification obtained from non-interacting fermions 
remains largely unknown. 

In this paper, wc fill this gap by considering 2D SCs 
with TRS {7^ = -1) and reflection symmetry (RS) 
{TZ^ = —1). (Note that T is anti- unitary while TZ is 
unitary.) We call SCs with these symmetry properties 
"DIII-I-R" SCs. First, we show that the topological 
classes of these DIII-I-R SCs are classified by an integer- 
valued (Z) topological invariant at the quadratic level. 
Such topologically non-trivial SCs are characterized by 
helical Majorana modes on the system's edge. The num- 
ber of helical Majorana edge states is protected by the 
symmetries at the non-interacting level. 

We then ask whether those helical gapless Majorana 
edge modes are stable against interactions while preserv- 
ing those relevant symmetries in question. As we show 
later in detail, the helical Majorana edge states in the 
case of Z = 8 topological invariant are unstable against 
(even weak) interactions - the gapless helical Majorana 
edge states become gapped while no symmetry is bro- 
ken in the bulk or the edges of the system. It is worth to 
stress that all the relevant symmetries are fully preserved 
while the gap opens in the edge. In other words, the edge 
of the putative Z = 8 system is qualitatively the same as 
Z = systems when symmetry-preserving interactions 
are included. A similar model, but with different set of 
symmetries, was studied independently in Refs. p4| . 

DIII+R SCs: Two-dimensional SCs with TRS 7^ = 
— 1 (symmetry class DIII) have a Z2 topological invariant 
as shown in the "periodic table" . To see this explicitly, it 
is helpful to study their edge theory. For a Z2 non-trivial 
DIII SCs in 2D, its generic edge theory is the ID helical 
Majorana fermions: 



(1) 



where ipa are left /right-moving edge Majorana fermion 
operators and a =t, 4- are spin index with the fol- 
lowing properties under time reversal transformation: 
T^^ipfT = ipi and T~^ipiT = — V'fj which satisfies = 
— 1. It is clear that the helical edge states above respect 
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the TRS. Moreover, this TRS protects this putative gap- 
less hehcal edge state since there is only one mass term 
imtlj^tp^ which breaks the TRS in the edge. However, 
when there are two such helical gapless modes labeled by 
a = 1,2, the putative gapless modes can be fully gapped 
by time reversal invariant terms im{1p^^p'^ + ipj^tp'^). Con- 
sequently, the topological classification is Z2. 

To possibly have a Z topological classification of Dili 
SCs in 2D, further discrete symmetry is needed. Since the 
helical edge states preserve not only TRS but also the RS 
TZ defined as Ti~^^^{x)TZ = ipil—x) and TZ~^ipi{x)TZ — 
-Vt(-2;) |25|, we consider 2D Dili SCs with RS, dubbed 
as the DIII-fR class. Note that TZ^ = — 1 as required 
for spin half-integer fermions. Now, we show that the 
DIII-I-R SCs in 2D are distinguished by an integer- valued 
(Z) topological invariant. To see this, we consider N 
copies (or flavors) of the helical gapless Majorana edge 

states given by X^^i / ^2; Va [iip^^dx'ip'f—ifpldxipl] , where 
a = 1 , . . . , are the "flavor" indices (which is preserved 
by both the time-reversal or reflection transformations) 
and Va are the Fermi velocities. We then write down the 
most general mass term 



(2) 



and check if it can preserve both TRS and RS. (Note 
that terms proportional to iil'ai'a =t or 4, are irrel- 
evant and only renormalize the Fermi velocities of edge 
Fermions.) Time-reversal symmetry requires Mab = 
while RS requires Mab = —M^a- Consequently, M = 
identically. In other words, it is impossible to gap the N 
pairs of helical gapless Majorana edge states by consid- 
ering the non-interacting fermion-bilinear terms for arbi- 
trary N . This class of SCs is then characterized by an 
integer- valued (Z) topological invariant (see below). 

Microscopic models: Armed with insights from the 
edge theories, we now introduce an explicit mean-field SC 
Hamiltonian which respects both TRS and RS and which 
has non-trivial gapless helical Majorana edge states: 

-ffiatt. = XI [ ~ ^^Lcjv + H.c] - fJ-J^^a^^'-^ 

{ij)a icr 
i 
i 

where are the fermion creation operators on site i 
and a i a-re spin indices. It is clear that the model 
above is invariant under RS IZ: Cx,y,a ^ {icf^)aa'C-x,y,a' 
and TRS T: c,^ ^ {i(^y)aa'Cra' This model de- 

scribes a SC with spin up Px +ipy pairing and spin down 
^Px + "ipy pairing, and can describe a thin film of '^He-B. 
On a cylinder with edges parallel to the x-direction, we 
(numerically) obtain the helical Majorana edge states, 
for i = A = 1 and for -4 < ^ < or < /x < -f4, 
say. These helical edge states are effectively described by 



the Hamiltonian (|l|). The TV- flavor generalization of this 
model is straightforward. 

Bulk topological invariant: There is a bulk topo- 
logical invariant which guarantees, when non-interacting 
and when there is a translation symmetry, the stability 
of the edge modes of DIII-I-R topological SCs for arbi- 
trary N. The construction is similar in spirit to the mir- 
ror Chern-number in 3D topological insulators protected 
by RS Q. With the periodic boundary condition, the 
quadratic bulk Hamiltonian can be Fourier transformed 



H 



(4) 



We then note at the reflection symmetric points kx = 
kx{= and tt), the Bloch Hamiltonian Hj, (ky) com- 
mutes with TZ, namely ['Hj.^{ky), J'^] = 0, where = 
diag(i(T^, — icT^). In other words, at kx, the quadratic 
Hamiltonian conserves S^. Combined with TRS, the 
quadratic Hamiltonian at kx = kx can be written as 



Hikx) = J2A+,k^<'-,-k)nkJky) 



(5) 



where the subscripts ± are eigenvalues of a^. Now, the 
Hamiltonian Til {ky) above is in the AIII class with chi- 
ral symmetry |14[|. 

Following Ref. gapped ID quadratic Hamiltonians 
in symmetry class AIII are distinguished by an integer 
topological invariant, the winding number v. Thus, at 
each reflection symmetric momentum kx = kx {kx = 
and tt), we can introduce an integer topological invariant 
("the reflection winding number"), v{kx), the winding 
number of "H^ {ky)- We thus have two integral topolog- 
ical invariants, v{kx = 0,7r). The non-zero value of the 
invariant, ^{kx) 7^ 0, guarantees the presence of |t'(fca:)| 
pairs of zero-energy Majorana states at kx, when an edge 
is introduced along x-direction. In particular, when the 
invariant is non-zero ^{kx) ^ at one of the reflection 
symmetric momenta {kx = 0, say) and it is zero at the 
other {kx = tt), this means there must be |i>(fc3,)| branches 
of non-chiral edge modes. In general, the difference of the 
reflection winding number 



v{kx = 0) - v{kx = tt) 



(6) 



tells us the number of non-chiral edge modes. 

Interaction effect: At the non-interacting level (and 
without disorder), the DIII-I-R SCs have Z topological 
classification as shown from both bulk and edge theories. 
Now, we consider the effect of symmetry-preserving in- 
teractions in the mean-field BdG Hamiltonian and check 
if the putative topological classification of Z is modified 
or not. Since gapless helical edge states are the hallmark 
of those topologically nontrivial SCs, we believe that it 
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would be sufficient to check if the gapless helical edge 
states are stable against interactions while requiring there 
is no symmetry breaking induced by interactions pTf . 

For the case of topological invariant N = 8 [more gen- 
erally N = (mod 8)], we try to identify certain in- 
teractions that can destabilize the gapless edge states 
while preserving the symmetries of the bulk and the edge. 
For simplicity, we consider the helical edge modes of free 
fermions with the same Fermi velocities Va = v 



actions 



(7) 



a=l 



which is invariant under the global SO (8) rota- 
tions among the left-moving or right-moving Majo- 
rana fermions. Now, we consider interactions allowed 
by TRS and RS. One naturally starts with the fol- 
lowing SO (8) symmetric interactions Hgn = -fffioc ^ 
g J dx(^ Sa=i ^'^'^i^i) I where g is the coupling constant. 
This theory is then the S0(8) Gross-Neveu (GN) model 
in (l-l-l)D, which is exactly solvable |2^. Especially, the 
interaction is marginally relevant for g > 0: for arbitrary 
small interaction strength g the ground state is gapped by 
spontaneously breaking the time reversal (or Z2 chiral) 
symmetry with the order parameter {iip'^'ip'^) ^ e~'^/^'"^\ 
There are twofold degenerate ground states at the edge. 
When g < 0, the interaction is marginally irrelevant 
and the ground state remains gapless. In other words, 
the SO (8) symmetric GN interactions cannot result in a 
unique gapped ground state in the edge. We need to look 
for some other channel of interactions to fulfill this. 

We follow the construction introduced in Rcf. [ pT[ . 
The non-interacting edge is described by the conformal 
field theory (CFT) of 8 free Majorana fermions, which 
is equivalent to the SO (8)1 Wess-Zumino-Witten model. 
The 8 fcrmion operators ■0^ (^ =t or J,) form the vector 
representation of SO (8). Moreover, the 16-dimensional 
spinor representation of SO (8) is reducible to two 8- 
dimensional irreducible ones formed by spinor operators 
77° and Xct, a, = 1, . . . , 8. The explicit forms of ry^ and 
are given by 



exp 



(8) 



where (j)'^ are boson fields obtained from bosonizing the 
system, ^iip'^'^ = e^""^" . The number of minus sign 

in the exponent of Eq. (||) is even for r/" but odd for 
Accidcntly, -0, rj and x ^-H form 8-dimension representa- 
tions of SO (8); they can actually be transformed into one 
another by the so-called triality symmetry of S0(8). It 
turns out that the spinor fields are useful to construct 
the interactions we desire. 

To fully gap the edge states without spontaneously 
breaking any symmetry, we consider the following inter- 



H,nt = - j dx ^(^EL/'^t'?! 



HO 



(9) 



which is S0(7)-invariant and leaves 77^ fixed. This S0(7)- 
symmetric interaction is also local in terms of the original 
fermions ip^. Indeed, a finite gap opens in the edge states 
while preserving the symmetries in question, as explicitly 
shown in Ref. |17| when B < and 2 A > B. To under- 
stand this, let us first look at the limit A ^ \B\ which is 
the SO (7) GN model plus free 77^ fermions (when B ^ 0). 
The chiral symmetry is broken by the SO (7) GN interac- 
tions with the order parameter M = {i X]I=i ^V^vl) 0: 
which generates a mass term iBMij^rj^ for 77® fermions. 
Now. the 77^ fermions can be mapped to the transverse 
field Ising model with field strength {h — he) oc zLBM, 
where he is the critical field strength in the transverse 
field Ising model which is equal to the Ising interaction 
strength. For h > he, its ground state is paramagnetic 
without any symmetry breaking; for h < he, the sys- 
tem spontaneously breaks the Ising symmetry resulting 
in twofold degenerate ground states. In other words, 
B > and i? < lie in two different phases. Since 
B = 2A > is the SO (8) GN model which sponta- 
neously break the chiral symmetry having twofold de- 
generate ground states, it is then clear that B < Q phase 
has a unique gapped ground state without breaking any 
symmetry. 

It is worthwhile to understand more heuristically why 
= 8 is special. The edge theory with = 8 is 
qualitatively equivalent to the two-leg ladder electron 
model at half-filling. Since there are two electrons per 
unit cell for the two-leg ladder at half-filling, having a 
fully gapped ground state without breaking any sym- 
metry is possible and expected To illustrate this, 
let us consider lattice Majorana fermions on Id chains 
described by the Hamiltonian H = uHi + WH2 
where Hi is quadratic in lattice real fermion operators: 
Hi = -(»/2)ELiE,SaA^"A^\i, where {X^,}^) - 25t^ 
and Sa = {1, 1, —1, —1, 1, 1, —1, —1}. It is clear that Hi 
is invariant under either TRS or RS, when we assume 



7^ : ^ (-l)''A!_+^A;+2 



-A°, 



-(-l)'^A^,-, (10) 
a= 1,2,5,6. (11) 



H2 is an interaction, and given by H2 = J2j W{X'j) where 
W{X'j) is a four fermion interaction composed of eight 
Majorana fermions: 

W{X'') = +X'X^X''X^ + A^A^A^A^ + A^A^A^A*^ (12) 
+ A^A^A^A^ + A^A^A^A^ + X^'X^'X'X' 

- A^A-^A'^A^ - AiA^A-'^A* + A^A^A^A^ + A^A^A^A^ 

- A^A^A^A^ - A^A^A^A^ - A^A^A^A^ - A^A^A^A^. 
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In the absence of the interaction term w = 0, the lat- 
tice model is gapless, whose continuum limit is given by 
the Hamiltonian (^. When we switch on lu ^ 0, the 
edge theory is gapped with unique ground state without 
breaking symmetries. This can be understood in the fol- 
lowing way: the interaction W can be written in terms 
of complex fermions 



{-\^ +iX^^)/2 = c. 



If 



2t' 



(A3 - jA4)/2 = ci^, 

(A^ + iA8)/2 = c2i, (13) 



as follows: 



WiX"^) = 165^,1 • Sj.2 + 2(71j.i - 1)2 + 2(71j.2 - 1)^ - 2, 

(14) 

where S'i_2 and n\ 2 is the spin and the fermion number 
operator for c\^2s (s =t / i)- With this interaction, the 
charge degrees of freedom will be frozen by Mott physics. 
The exchange interaction S\- S2 (the "rung-exchange" 
interaction) realizes the rung-single phase which has a 
unique ground state without breaking any symmetry. 

Discussion: We have discussed interaction effects on 
topological SCs in 2D protected by TRS and RS. This 
is a non-trivial example in 2D where topological classi- 
fication of ground states are dramatically altered by the 
interaction effect. In a separate paper, we plan on a 
more systematic study on topological insulators and SCs 
protected by RS in addition to other possible discrete 
symmetries |3^, |3^ . 

We close with discussion on disorder effects: In the 
ten-fold classification of topological insulators and SCs, 
it has been proved useful to consider the boundary (edge, 
surface, etc.) Anderson localization problem: For a topo- 
logical bulk, one should find a boundary mode which 
is completely immune to disorder. In turn, once one 
finds such "Anderson delocalization" at the boundary, it 
means there is a topologically non-trivial bulk. Not only 
this bulk-boundary correspondence can be used to find 
and classify bulk topological phases in the absence of dis- 
order, it immediately tells us such topological phases are 
stable against disorder. For topological phases protected 
by a set of spatial symmetries, stability against disor- 
der is, in general, not trivial, since spatial inhomogeneity 
does not respect the spatial symmetries. One can still 
consider, however, situations where the spatial symme- 
tries are preserved on average. The effects of disorder in 
TV-channel quantum wires in symmetry class Dili, which 
are from our point of view the edge theory of the DIII+R 
topological SC with the topological integer (= N), have 
been studies 32 . It is known that there is an even-odd ef- 



fect in N: the mean conductance decreases algebraically 
as L~^/2 with the length of the wire L for odd N, whereas 
it decays exponentially with L for even N. Correspond- 
ingly, for symmetry class Dili with N odd the density of 
states the Dyson singularity. This implies that disorder 
simply reduces the Z topological classification of DIII+R 



SCs to the Z2 classification, which is the same as the 
topological classification of 2D Dili SCs in the periodic 
table. In a separate paper we will report the other 
cases where the Z topological classification is reduces to 
the Z2 classification, which is not related to the existing 
topological class in the periodic table. 

Acknowledgements: Wc would like to thank Eun- 
Ah Kim, Steve Kivelson, Dung-Hai Lee, and Shou-Cheng 
Zhang for helpful discussion, and in particular Xiao- 
Liang Qi for sharing his results with us prior to its arXiv 
submission. This work is supported in part by NSF Grant 
DMR-0904264. 



[1] 
[2] 
[3] 
[4] 

[6: 

[7 

[8 

[9: 
[10: 

[11 

[12 



[is: 

[14 

[is: 

[16: 
[17: 

[is: 

[19 

[20 

[21 

[22: 

[23' 
[24: 



The Quantum Hall Effect, edited by R. E. Prange and S. 

M. Girvin (Springer, New York, 1987). 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011) . 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005); C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 
226801 (2005). 

B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 

106802 (2006). 

R. Roy, Phys. Rev. B 79, 195321 (2009). 

J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) 

(2007). 

L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 

R. Roy, Phys. Rev. B 79, 195322 (2009). 
L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
X.-L. Qi and T. Hughes and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008). 

See, Y. Wada et al, Phys. Rev. B 78, 214516 (2008), S. 
Murakawa et al, Phys. Rev. Lett. 103, 155301 (2009), S. 
Murakawa et al, J. Phys. Soc. Jpn. 80, 013602 (2011), 
and references therein. 

Here, a "non-interacting" SC means it is described by its 
mean-field Hamiltonian even though the SC comes from 
some sort of interactions. 

A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 

wig, Phys. Rev. B 78, 195125 (2008). 

S. Ryu, A. Schnyder, A. Furusaki and A. W. W. Ludwig, 

New J. Phys. 12, 065010 (2010). 

A. Kitaev, AIP Conf. Proc. 1134, 22 (2009). 

L. Fidkowski and A. Kitaev, Phys. Rev. B 81, 134509 

(2010) . 

L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 

(2011) . 

A. M. Turner, F. Pollmann, and E. Berg, Phys. Rev. B 
83, 075102 (2011). 

S. Ryu, J. E. Moore, and A. W. W. Ludwig, Phys. Rev. 
B 85, 045104 (2012). 

Z. Wang, X.-L. Qi, and S.-C. Zhang, Phys. Rev. B 84, 
014527 (2011). 

X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, 

arXiv:1106.4772 (2011). 

X.-L. Qi, arXiv: 1202.3983 (2012). 

S. Ryu and S.-C. Zhang, arXiv:1202.4484 (2012). 



5 



[25] Note that the reflection TZ is given by ia^ in the bulk. 
When projected on tlie edge, TZ is given by ia^ , as re- 
quired by the reality nature of Majorana fermions on the 
edge and the condition Ti? = — 1. 

[26] J. C. Y. Teo, L. Fu, and C. L. Kane, Phys. Rev. B 78, 
045426 (2008). 

[27] To explicitly show that the topological classification in 
the interacting DIII+R SCs is Zg, constructing an adi- 
abatic path connecting the N = 8 model to the = 
model without closing the bulk gap should be needed, as 
what Fidkowski and Kitaev did in [|I^ for the (l-fl)D 
BDI SCs, even though we believe that the classification 



argued from the edge consideration is sufficient. 
[28] R. Shankar, Phys. Rev. Lett. 46, 379 (f98f). 
[29] H.-H. Lin, L. Balents, and M. P. A. Fisher, Phys. Rev. B 

58, 1794 (1998). 
[30] H. Yao and S. Ryu, in preparation. 

[31] A systematic study of ID edge theories of a 2D bulk was 

also given recently in D. Bernard, E. -A. Kim, and A. 

LeClair, arXiv: 1202.5040. 
[32] See, for example, P. Brouwer, A. Furusaki, C. Mudry, and 

S. Ryu, BUTSURI 60, 935 (2005) (in Japanese) (English 

version: arXiv:cond-mat /05 1 1622) . 



